Magnetism of CuX 2 frustrated chains (X = F, CI, Br): the role of covalency 
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Periodic and cluster density-functional theory (DFT) calculations, including DFT+£7 and hybrid 
functionals, are applied to study magnetostructural correlations in spin-1/2 frustrated chain 
compounds CuX2: C11CI2, CuBr2, and a fictitious chain structure of CuF2. The nearest-neighbor 
and second- neighbor exchange integrals, J\ and J2, are evaluated as a function of the Cu-X-Cu 
bridging angle in the physically relevant range 80-110°. In the ionic CuF2, Ji is ferromagnetic for 
< 100°. For larger angles, the antiferromagnetic superexchange contribution becomes dominant, 
in accord with the Goodenough-Kanamori-Anderson rules. However, both C11CI2 and CuBr 2 
feature ferromagnetic Ji in the whole angular range studied. This surprising behavior is ascribed 
to the increased covalency in the CI and Br compounds, which amplifies the contribution from 
Hund's exchange on the ligand atoms and renders J\ ferromagnetic. At the same time, the larger 
spatial extent of X orbitals enhances the antiferromagnetic J2 , which is realized via the long-range 
Cu-X-X-Cu paths. Both, periodic and cluster approaches supply a consistent description of the 
magnetic behavior which is in good agreement with the experimental data for C11O2 and CuBr2. 
Thus, owing to their simplicity, cluster calculations have excellent potential to study magnetic 
correlations in more involved spin lattices and facilitate application of quantum-chemical methods. 
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I. INTRODUCTION 

Copper compounds have been extensively studied as 
spin- ^ quantum magnets, material prototypes of quan- 
tum spin models. While local properties of these com- 
pounds are usually similar and involve nearly isotropic 
Heisenberg spins, the variability of the magnetic behavior 
stems from the unique structural diversity. Depending on 
the particular arrangement of the magnetic Cu 2+ atoms 
and their ligands in the crystal structure, different spin 
lattices can be formed. Presently, experimental examples 
for many of simple lattice geometries, including the uni- 
form chain, 1 ' 2 square lattice, 3 ' 4 Shastry- Sutherland lat- 
tice of orthogonal spin dimers, 5 ' 6 and kagome lattice, 7 
are available and actively studied. Some of the copper 
compounds feature more complex spin lattices 8 10 that 
have not been anticipated in theoretical studies, yet trig- 
ger the theoretical research 11 ' 12 once relevant material 
prototypes are available. 

Owing to the competition between ferromagnetic (FM) 
and antiferromagnetic (AFM) contributions to the ex- 
change couplings, compounds of particular interest are 
those with M-X-M bridging angles close to 90°, with M 
being a transition metal and X being a ligand. Such 
geometries are realized in the quasi- ID cuprates featur- 
ing chains of edge-sharing CUO4 plaquettes, which repre- 
sent a simple example of low-dimensional spin-1/2 mag- 
netic materials. Independent of the sign of the nearest- 
neighbor (NN) coupling Ji, its competition with the 
sizeable AFM next-nearest-neighbor (NNN) coupling J2 
leads to magnetic frustration. Depending on the ratio 
J2/J11 such compounds exhibit exotic magnetic behavior 
like helical order, 13 spin-Peierls transition 14 or quantum 
critical behavior. 15 The difficulties in the microscopic de- 



scription of such compounds originate from ambiguities 16 
in the experimental estimates of the ratio J2/J1, lead- 
ing to controversial modeling of the magnetic structure. 17 
Thus, the combination of different sets of experimental 
data with a careful theoretical analysis of the individual 
exchange pathways is of crucial importance for obtaining 
a precise microscopic magnetic model. 

However, the search for new quantum magnets, as well 
as the work on existing materials, require not only the 
ability to estimate the couplings but also a solid under- 
standing of the nexus between crystallographic features 
of the material and ensuing magnetic couplings. The 
Goodenough-Kanamori-Anderson (GKA) 18 rules are a 
generic and well-established paradigm that prescribes 
FM couplings for bridging angles close to 90° and AFM 
couplings else, where the bridging angle refers to the M- 
X-M pathway. In Cu 2+ oxides, generally the GKA rules 
successfully explain the crossover between the FM and 
AFM interactions for Cu-O-Cu angles close to 90°. The 
boundary between the FM and AFM regimes is usually 
within the range of 95 — 98°, 19 but may considerably be 
altered by side groups and distortions. 20 ' 21 

In addition to Cu 2+ oxides, the systems of interest 
include copper halides, 22 ' 23 carbodiimides, 24 and other 
compound families. Although microscopic arguments be- 
hind the GKA rules should be also applicable to these 
non-oxide materials, the critical angles separating the 
FM and AFM regimes, as well as the role of the ligand 
in general, are still little explored. Moreover, the low 
number of experimentally studied compounds impedes a 
comprehensive experimental analysis available for oxides. 

More and more, density functional theory (DFT) 
electronic structure calculations complement experimen- 
tal studies and deliver accurate estimates of magnetic 
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couplings. 25 30 They are especially well suited for the 
study of magnetostructural correlations, as both real 
and fictitious crystal structures can be considered in a 
calculation. However, in a periodic structure the ef- 
fect of a single geometrical parameter is often difficult 
to elucidate, because different geometrical features are 
intertwined and evolve simultaneously upon the varia- 
tion of an atomic position. Geometrical effects on the 
local magnetic coupling are better discerned in cluster 
models that represent a small group of magnetic atoms 
and, ideally, a single exchange pathway. Additional ad- 
vantages of cluster models, owing to their low number 
of correlated atoms, are lower computational costs and, 
most important, their potential for the application of 
parameter-free wavefunction-based computational meth- 
ods, i.e. in a strict sense ab initio calculations. By con- 
trast, presently available band- structure methods for cal- 
culating strongly correlated compounds rely on empirical 
parameters and corrections where their choice is in gen- 
eral not unambiguous. 31,32 

There have been several attempts to describe the lo- 
cal properties of solids with clusters especially in com- 
bination with ab initio quantum-chemical methods. 33 37 
However, the construction of clusters is far from being 
trivial. On one side, to make the calculations compu- 
tationally feasible, the number of quantum mechanically 
treated atoms has to be kept as small as possible. On 
the other side, accurate results require that these atoms 
experience the "true" crystal potential. Usually, this is 
achieved by embedding the cluster into a cloud of point 
charges 36,38 and so called total ion potentials. 39,40 But 
even for involved embeddings it was demonstrated, that 
the choice of the cluster may have significant effects on 
the results of the calculations and, thus, size-convergence 
has to be checked thoroughly. 40 

Here, we study the effect of geometrical parameters 
on the magnetic exchange in Cu 2+ halides. The family 
of halogen atoms spans a wide range of electronegativi- 
ties, from the ultimately electronegative fluorine, form- 
ing strongly ionic Cu-F bonds, to chlorine and bromine 
that produce largely covalent compounds with Cu 2+ . 41 
Presently, we do not consider iodine because no Cu 2+ io- 
dides have been reported. In our modeling, we use the 
simplest possible periodic crystal structure of a CuX2 
chain that enables the variation of the Cu-X-Cu bridg- 
ing angle in a broad range. We further perform a com- 
parative analysis for clusters and additionally consider 
the problem of long-range couplings. The evaluation of 
such couplings requires larger clusters, thus posing a dif- 
ficulty for the cluster approach. The observed trends 
for the magnetic exchange as a function of the bridging 
angle are analyzed from the microscopic viewpoint, and 
reveal the crucial role of covalency that underlies salient 
differences between the ionic Cu 2+ fluorides and largely 
covalent chlorides and bromides. 

On the experimental side, the compounds and crys- 
tal structures under consideration are relevant to the 
CuCl2 and CuBr2 materials that show interesting ex- 



amples of frustrated Heisenberg chains. 42,43 At low tem- 
peratures, these halides form helical magnetic structures 
and demonstrate improper ferroelectricity along with the 
strong magnetoelectric coupling. 44,45 

The paper is organized as follows. In section II, the 
applied theoretical methods are presented. In the third 
section, the crystal structures of the CuX2 compounds 
are described and compared. In section IV, the results of 
periodic and cluster calculations are discussed and com- 
pared. Finally, the discussion, summary, and a short 
outlook are given in section V. 



II. METHODS 

The electronic structures of clusters and periodic sys- 
tems were calculated with the full-potential local-orbital 
code fplo9. 00-34. 46 For the scalar-relativistic calcu- 
lations within the local density approximation (LDA), 
the Perdew-Wang parameterization 47 of the exchange- 
correlation potential was used together with a well con- 
verged mesh of up to 12x12x12 k-points for the periodic 
models. 

The effects of strong electronic correlations were con- 
sidered by mapping the LDA bands onto an effective 
tight-binding (TB) model. The transfer integrals t{ of 
the TB-model are evaluated as nondiagonal elements be- 
tween Wannier functions (WFs). For the clusters, the 
transfer integral corresponds to half of the energy differ- 
ence of the magnetic orbitals. 48 These transfer integrals 
U are further introduced into the half-filled single-band 
Hubbard model H = Htb + ^eff n it n ii that is even- 
tually reduced to the Heisenberg model for low-energy 
excitations, 

H = Y,JiA-Sj (i) 

(ij) 

The reduction is well-justified in the strongly correlated 
limit ti <C E/eff , where f7 e ff is the effective on-site Coulomb 
repulsion, which exceeds U by at least an order of mag- 
nitude (see Table I). This procedure yields AFM contri- 
butions to the exchange evaluated as jf^ FM = 4tf/U e f[. 

Alternatively, the full exchange couplings J^, compris- 
ing FM and AFM contributions, can be derived from 
total energies of collinear magnetic arrangements eval- 
uated in spin-polarized supercell calculations 49 within 
the mean-field density functional theory (DFT)+Z7 for- 
malism. We use a local spin-density approximation 
(LSDA)+/7 scheme in combination with a unit cell 
quadrupled along the b axis and a /c-mesh of 64 
points. The on-site repulsion and exchange amount to 
Ud = 7±0.5 eV and Jd = 1 eV, respectively. The same Ud 
value is chosen for all CuX2 (X = F, CI, Br) compounds to 
facilitate a comparison of the magnetic behavior. In sec- 
tion IV E, however, it will be shown that Ud has in fact no 
qualitative effect on the magnetic couplings of the CuX2 
(X = F, CI, Br) compounds. We applied the around mean 



field (AMF) as well as the fully localized limit (FLL) dou- 
ble counting corrections where both types where found to 
supply similar results. Thus, following the earlier stud- 
ies of Cu 2+ compounds, 25,26,43 the presented results are 
obtained within the AMF scheme. 

For the clusters we used, in addition to the LSDA+/7 
method, the B3LYP hybrid functional 50 with a 6-311G 
basis set. The B3LYP calculations were performed within 
the GAUSSIAN09 code. 51 The free parameter a, indicating 
the admixture of exact exchange, was varied in the range 
between 0.15 and 0.25 to investigate its influence on the 
calculated exchange couplings. 

III. CRYSTAL STRUCTURES 
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FIG. 1. (Color online) Edge-sharing CuX4-plaquettes form- 
ing the magnetic chains in the CuX2 compounds. The chains, 
running along [010] are flat and lie in the ab plane. The 
stacking of the planes is accompanied by a shift to match the 
monoclinic angle. 52 The arrows indicate the nearest-neighbor 
and next-nearest-neighbor interaction pathways, and de- 
notes the Cu-X-Cu bridging angle. 



The copper CuX2 dihalides feature isolated chains 
of edge-sharing CUX4 plaquettes. 52 The chains of this 
type are the central building block of many well-studied 
cuprates such as CuGe03 (Ref. 14), Li2ZrCu04 (Ref. 15), 
and Li 2 Cu02 CuX 2 halides are charge neutral, which 
makes them especially well suited for the modeling within 
the cluster approach. 

CuBr2 crystallizes in the monoclinic space group 
C2/m with a = 14.728 A, b = 5.698 A and c = 8.067 A, and 
f3 = 115.15° at room temperature. 53 The planar chains 
of edge-sharing CuBr4 plaquettes run along the 6-axis 
(Fig. 1). The Cu-Br-Cu bridging angle amounts to 
92.0°, the Cu-Br distance is 2.41 A, while the distances 
between the neighboring chains amount to dy=3.82A 
and d± = 3.15 A in the direction parallel to c and per- 
pendicular to the plaquette plane, respectively. 

CuCl2 is isostructural to CuBr 2 with the Cu-Cl dis- 
tance of 2.26 A and Z(Cu-Cl-Cu) = 93.6°. 54 The inter- 
chain separations amount to d\\ = 3.73 A and d± = 2.96 A 
along the c and a directions, respectively. 

CuF 2 features a two-dimensional distorted version of 
the rutile structure, with corner-sharing CuF 4 plaque- 
ttes forming a buckled square lattice. 55 This atomic ar- 
rangement is very different from the chain structures of 
CuCl2 and CuBr2. For the sake of comparison with other 
Cu 2+ halides, we constructed a fictitious one-dimensional 
structure of CuF2. The Cu-F distance of 1.91 A was cho- 
sen to match the respective average bond length in the 
real CuF 2 compound. The corresponding bridging angle, 
yielding a minimum in total energy, was determined to 
be 102°. 56 Although this crystal structure remains hy- 
pothetical, it is likely metastable and could be formed 
in CuF 2 under a strong tensile strain on an appropriate 
substrate. 



IV. RESULTS 
A. Band structure calculations 

First, we consider magnetic couplings in the experi- 
mental crystal structures of CuCl2 and CuBr2, as well as 



in the relaxed structure of chain-like CuF2- The DFT cal- 
culations of the band structure and the density of states 
(DOS) of CuX 2 (X = F, CI, Br) compounds within the 
LDA yield a valence band width of 6-8eV, 52 in agree- 
ment with the experimental photoelectron spectra. 41 The 
valence band complex becomes slightly narrower upon an 
increase in the ligand size, because the lower electronega- 
tivity of CI and Br brings the respective p states closer to 
the Cu 3d states, thus enhancing the hybridization and 
reducing the energy separation between the Cu and lig- 
and orbitals. All the band structures feature a separated 
band crossing the Fermi level (Fig. 2). In the local-orbital 
representation visualized by WFs (Fig. 4), this band is 
formed by the antibonding <j*-combination of Cu 3d x 2_ y 2 
and X p orbitals. 57 The isolated half-filled band suffices 
for describing the magnetic properties and the low-lying 
magnetic excitations via the transfer integrals t{ which 
are subsequently introduced into a Hubbard model. Lig- 
and valence p-orbital contributions to the magnetic or- 
bital, denoted as j3 in Table I, illustrate the increase in 
the metal-ligand hybridization from F to Br. 

The dispersion calculated with the WF-based one- 
band TB model for CuBr2 is also shown in Fig. 2, and 
the leading transfer integrals together with the AFM con- 
tributions J 2 AFM are given in Table I. The evaluation of 
jAFM re q U i res the value of J7 e ff, which is not known pre- 
cisely. Here, we estimate f7 e ff by comparing the trans- 
fer integral £2 obtained from the TB analysis with the 
exchange coupling J2 from the LSDA+Z7 calculations. 
While short-range couplings may involve large FM contri- 
butions, the long-range coupling J 2 should be primarily 
AFM. Therefore, J AFM = J2 in a good approximation, 
and U e ?E = 4£ 2 / J 2- This way, we find J7 e ff = 6 eV for X = 
F, 4 eV for CI, and 3 eV for Br. The reduction in Z7 e ff re- 
flects the general trend of the enhanced Cu-X hybridiza- 
tion and covalency, because the f7 e ff value pertains to the 
screened Coulomb repulsion in the mixed Cu-X band. 
The enhanced hybridization leads to a stronger screen- 
ing, larger spatial extension and, thus, to the lower f7 e ff 
values. 

The estimates in Table I reveal two major differences 



4 



TABLE I. Results for the experimental (X = CI, Br) and hypothetical (X = F) structures of CUX2: the bridging angle 
the ligand contribution to the magnetic orbital /3, transfer integrals U, AFM contributions to the exchange J 2 AFM = 4tf/U e s, 
total exchange integrals Ji from LSDA+£7 calculations with Ud — 7 =b 0.5 eV, and the effective on-site Coulomb repulsion U e s 
obtained by equilibration of At^/Ues (see text for details). 
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Finally, we briefly compare our DFT-based estimates 
of Ji with the experiment. Because the chain-like poly- 
morph of CuF2 has not been prepared experimentally, 
no comparison can be performed. The microscopic anal- 
ysis of CuCl2 presented in Ref. 43 shows reasonable 
agreement between the experimental (Ji = —7.8 meV, 
J 2 = 11.6 meV) and calculated (Ji = -12.9 ± 0.9 meV, 
J 2 = 13.4 ±2.2 meV) values. The same is true for 
CuBr 2 , where we evaluated the intrachain couplings as 
Ji = -8.8 ± 0.4 meV, J 2 = 22.2 ± 3.4 meV which 
compare well with recently published experimental data 
Ji = -11.0 ± 1.6 meV, J 2 = 31.0 meV. 58 Moreover, our 
calculations reveal significantly lower deviations from ex- 
periment than those supplied in Ref. 58. 



FIG. 2. (Color online) Comparison of the calculated LDA 
band structure of CuBr 2 and the band derived from a fit us- 
ing an effective one-band tight-binding model based on Cu- 
centered Wannier functions (WF TB). The right plot shows 
the total density of states (DOS) together with the partial 
DOS of Cu(3cf) and Br(4p) states. The Fermi level is at zero 
energy. Notation of fc-points: T = (000), X= (-00), S= 
(*f 0), Y= (Of 0), YZ= (Off), XYZ= (I f f ), Z= a ((H)f ). 



between the ionic CuF 2 and more covalent CuCl 2 and 
CuBr 2 compounds. First, the nearest-neighbor (NN) 
coupling Ji is AFM in the fluoride, while FM in the 
chloride and bromide. Second, the AFM next-nearest- 
neighbor (NNN) coupling J 2 is enhanced upon increasing 
the covalency of the Cu-X bonds. In CuF 2 , this coupling 
is weak (J 2 <C Ji), whereas in the chloride and bromide 
J2 > |Ji|- The NNN coupling is amplified by the larger 
ligand size and the increased covalency. This coupling in- 
volves the long-range Cu-X-X-Cu pathway and requires 
a strong overlap between the ligand orbitals, which is 
possible for X = CI and especially Br, while remaining 
weak for the smaller fluoride anion. The changes in the 
NN coupling seem to be well described by the GKA rules. 
Considering the trends for copper oxides, 19 one expects 
FM Ji for 6 close to 90°, as in CuCl 2 and CuBr 2 , and 
AFM Ji for 6 > 98°, as in the chain-like structure of 
CuF 2 . Nevertheless, the covalency is also paramount for 
the sign of Ji, as shown by the magnetostructural corre- 
lations presented below (Sec. IV B). 



Puzzled by the origin of the discrepancy between our 
values for Ji and J 2 and the published calculational re- 
sults for CuBr 2 , 58 we repeated the DFT±L^ calculations 
for CuBr 2 as well as CuCl 2 with the code VASP 59 and 
the same computational parameters as used in Ref. 58. 
For the parameters Ud and J^, we adopted 8eV and 
1 eV, respectively, which corresponds to the effective 
U= U d -J d = 7eV in Refs. 58. For the GGA+Z7 cal- 
culations, we used again a unit cell quadrupled along the 
b axis and the /c-mesh of 64 points. The resulting Ji and 
J 2 values generally agree with the published values, 42,58 
except for Ji in CuBr 2 , for which we obtain only half of 
the value provided in Ref. 58. The agreement with the 
experimental data can be improved by increasing the Ud 
value. In particular, U d = 12 eV yields Ji = — 95 K and 
J 2 = 113K for CuCl 2 and J 1 = -124 K and J 2 = 357K 
for CuBr 2 , very close to the experimental estimates. 43,58 
This Ud value is significantly higher than the U d = 7eV 
we used in our fplo9. 00-34 calculations. 60 There are 
basically two reasons for the large difference: The first 
reason are the different basis sets of fplo9. 00-34 and 
VASP, implementing local orbitals and projected aug- 
mented waves, 61 respectively, which crucially affect the 
local quantity U d - Second, we used an around mean field 
double counting correction (DCC) while a fully localized 
limit DCC, which is always used in VASP, requires larger 
Ud values. 62 
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FIG. 3. (Color online) J\ of CuCb as function of the bridging 
angles where different structural parameters are fixed: i) the 
Cu-Cu distance is varied, while the X position is subsequently 
optimized to yield the equilibrium Cu-X distance and 0; ii) 
the Cu-Cu distance is fixed, while the Cu-X distance is varied; 
and iii) the Cu-X distance is fixed, while the Cu-Cu distance 
is varied. The dashed vertical line indicates the experimental 
bridging angle. 



B. Variation of the bridging angle 

To establish magnetostructural correlations in CuX2 
halides, we systematically vary the bridging angle and 
evaluate the NN coupling J\. Since the Cu-Cu distance 
and two Cu-X distances form a triangle with being one 
of its angles, the change in alters either the Cu-Cu 
distance, or the Cu-X distance, or both. We compared 
different flavors of varying 0: 6S i) the Cu-Cu distance is 
varied, while the X position is subsequently optimized to 
yield the equilibrium Cu-X distance and #; ii) the Cu- 
Cu distance is fixed, while the Cu-X distance is varied; 
and iii) the Cu-X distance is fixed, while the Cu-Cu dis- 
tance is varied. For all three cases, we evaluated J\ as a 
function of the Cu-X-Cu angle. Fig. 3 shows on the ex- 
ample of CuCl2 that despite minor numerical differences, 
all three methods conform well to each other. Addition- 
ally, we studied the influence of Ud by varying it in the 
wide range of 4-9 eV. This causes a shift of the curves 
along the vertical axis, but the qualitative behavior of J\ 
versus the Cu-X-Cu angle is retained. 52 

Remarkably, J\ reaches its minimum absolute value 
at around = 100° and becomes strongly FM at large 
bridging angles (Fig. 3). This result is robust with re- 
spect to the particular procedure of varying 0. To better 
understand the microscopic origin of this peculiar be- 
havior, we performed similar calculations for CuF2 and 
CuBr 2 . As different procedures of varying arrive at sim- 
ilar results, we fixed the Cu-X distance for each ligand 
and achieved different values by adjusting the Cu-Cu 
distance, only. 



Similar to our results for the fixed geometries (Table I) , 
magnetostructural correlations for J\ (Fig. 4) reveal a 
large difference between the ionic CuF 2 and covalent 
CuCl2 and CuBr2. In CuF2, J\ follows the anticipated 
behavior with the FM-to-AFM crossover at ~ 100°. 
However, the covalent compounds always show FM Ji, 
with a maximum (i.e., the minimum in the absolute 
value) at = 100 — 105° and the enhanced FM charac- 
ter at even larger bridging angles. This trend persists up 
to at least = 120° (Fig. 3). 

The effect of strongly FM J\ in CuCl2 and CuBr2 can 
be explained by considering individual contributions to 
the exchange. The AFM contribution J AFM arises from 
the electron hopping between the Cu sites. The hopping 
probability measured by the transfer integral t\ critically 
depends on the Cu-X-Cu bridging angle. In a simple 
ionic picture, the transfer is maximal at = 180° (singly 
bridged) and approaches zero at = 90°, thus provid- 
ing the microscopic reasoning behind the GKA rules. 
This anticipated trend is indeed shown by CuF2, where 
jAfm _ 4tf /f/ e ff increases above = 90° and underlies 
the increase in J\. However, the covalent CuX2 halides 
show qualitatively different behavior with the very low 
(and decreasing) t\ and J AFM up to at least — 110°. 
This result implies that the large contribution of the lig- 
and states in a covalent compound has also a strong in- 
fluence on the Cu-X-Cu hopping process and alters the 
anticipated trend for the AFM exchange. 

The FM contribution J FM can be evaluated as 
J\ — J AFM , where we use J\ from the LSDA+/7 calcu- 
lation and J X AFM = Atl/Ues from the TB analysis. Mi- 
croscopically, J FM originates from the Hund's coupling 
on the ligand site 64 and/or from the FM coupling be- 
tween the Cu 3d and ligand p states. 65 Regarding the 
former mechanism, 64 a simple model expression reads as 
jFM _ _^4j^ w here f3 is the ligand's contribution to 
the Cu-centered magnetic orbital, and Jh is the (effec- 
tive) Hund's coupling on the ligand. Even though this 
expression is derived for = 90°, our data obtained for 
different values are well understood in terms of the vari- 
able (3 (see bottom panels of Fig. 4). The increase in the 
bridging angle leads to larger /?, thus enhancing J FM . 
Since f3 enters Jf M as /3 4 , its effect should be dominant 
over any other contributions, such as slight variations of 
Jh- The increase in j3 also explains the increasing FM 
contribution at low (Fig. 4). 

In contrast to the covalent chloride and bromide, the 
ionic CuF2 shows only a minor FM contribution owing 
to the very low j3. We also tried to artificially enhance /3 
by reducing the Cu-F bonding distance down to 1.60 A. 
For bridging angles larger than 100° the AFM coupling 
becomes twice as large as for the Cu-F distance of 1.91 A 
and for angles smaller than 80° the model compound be- 
comes also AFM. The FM coupling strength about 90° 
is almost unaffected. This indicates the robust ionic na- 
ture of Cu-F bonds. The reduction in the Cu-F distance 
increases the electron transfer without changing the hy- 
bridization, hence J AFM is increased, while Jf M remains 
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weak. 



C. Cluster models 

In a periodic calculation, the variation of structural pa- 
rameters, such as bond lengths and angles, is generally 
challenging: the high symmetry couples the structural 
parameters to each other. As a result, changing a sin- 
gle parameter is often impossible without affecting the 
other parameters. The cluster models are more flexible 
and may allow for an independent variation of individ- 
ual bond lengths and angles. This property renders the 
clusters as an excellent playground to study the magne- 
tostructural correlations. 

Before discussing the intrachain couplings using a com- 
bination of periodic and cluster models, we first want 
to demonstrate how cluster models for the three Cu di- 
halide compounds are constructed. Since the chains are 
spatially well- separated from each other, we can consider 
segments of a chain, with the terminal Li atoms keep- 
ing the electroneutrality (Fig. 5). No additional point 
charges are required, so that the clusters are kept as sim- 
ple as possible. 

First, the effect of the chain length on Ji, J 2 and the 
ratio —J2/J1 is investigated (Fig. 6). For all three com- 
pounds, small clusters, such as dimers or t rimers, are 
insufficient for describing the magnetic properties. The 
convergence with respect to the cluster size is different for 
different compounds (e.g., the ionic CuF 2 demonstrates 
the slowest size convergence). To ensure a meaningful 
comparison with the periodic model or the experimental 
data, the convergence with respect to the cluster size has 
to be carefully checked. 

On the other hand, a large number of correlated centers 
requires a large number of spin configurations to estimate 
exchange couplings. While larger clusters are still feasible 
for DFT, they may pose a problem for advanced ab initio 
quantum-chemical methods. Therefore, we attempted to 
reduce the number of correlated Cu 2+ ions by substitut- 
ing them by formally nonmagnetic Mg 2+ and Zn 2+ ions 
(Fig. 5). Even with this minimum number of correlated 
centers, deviations below 10% to the size-converged Cus 
octamer cluster are obtained for the Cu-Br (Fig. 7) and 
also for the Cu-Cl clusters. In case of Cu-F, where con- 
vergence is reached at larger cluster size, at least four 
correlated centers are required to reduce the deviations 
down to that level. 

Similar results, as for the J's, concerning size con- 
vergence and substitutions are obtained for the NN and 
NNN transfers, t\ and £2, calculated in LDA. These re- 
sults show that the simple clusters suffice for describ- 
ing the intrachain physics of these compounds and that 
the problem of appropriately embedding the clusters may 
be at least partially bypassed by increasing the cluster 
size and substituting part of the correlated centers with 
weakly correlated ions. 



D. Cluster versus periodic models 

In the following, both cluster and periodic models will 
be used for calculating J 2 and the —J2/J1 ratio, as well 
as the transfer integrals t{ of the Cu dihalides. The com- 
parison of periodic and cluster models for a broad range 
of bridging angles allows to exclude an accidental agree- 
ment between both models, which can be realized in a 
specific geometry by appropriately choosing the chain 
length, substitutions, and the termination of the cluster. 
However, when the cluster is prepared in such a way, the 
good agreement with the periodic model would be lost 
by varying the geometrical parameters. 

The exchange integrals as well as the —J2/J1 ratio ver- 
sus the bridging angles are depicted in Figs. 8 and 9 for 
CuBr2 and CuCl2, respectively. A comparison of the 
nearest-neighbor transfer integral ti of CuBr 2 , calculated 
with cluster and periodic models, is shown in Fig. 10. 
The clusters can reproduce the results of band structure 
calculations over the whole range of bridging angles, thus 
justifying the construction of the clusters. In the —J2/J1 
ratio, which governs the magnetic ground state, the devi- 
ations between the cluster and periodic models are com- 
pensated to a large degree. 

In CuF2, the deviations between J\ and J 2 obtained 
in the cluster and periodic models, respectively, are also 
compensated in the ratio — J2/J1, except for the smallest 
bridging angles. 52 The singularity in — J2/J1 at about 
100° arises from the crossover between the FM and the 
AFM Jx. 

These results show that well-controlled cluster mod- 
els are capable of describing local properties of ionic as 
well as strongly covalent solids, whereas the good agree- 
ment with band structure calculations is not accidental or 
artificial. Finally, the results demonstrate that superex- 
change and magnetic coupling in insulators are relatively 
short-range effects even for strongly covalent compounds. 

E. LSDA+?7 vs. hybrid functionals 

A common problem of DFT-based approaches applied 
to strongly correlated electrons is the ambiguous choice of 
empirical parameters and corrections that are required to 
mimic many-body effects, e.g., in the mean- field DFT+U 
approach. Hybrid functionals represent an alternative, 
although still empirical, way of simulating the effect of 
strong electron correlations within DFT. In this way, the 
non-local exact exchange is mixed with the local LDA 
or GGA exchange, while the mixing parameter a is typ- 
ically the only free parameter. In contrast to DFT+C7, 
hybrid functionals are more robust with respect to the ad- 
justable parameters, and the constant value of a — 0.20 
or a = 0.25 can be used in a rather general fashion. 
Additionally, the exact exchange correction is generally 
applied to all orbit als while in DFT+/7 the corrections 
are applied to a certain set of orbitals which are assumed 
to be the strongly correlated ones. 
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FIG. 4. (Color online) Magnetostructural correlations for the CuX2 halides (the Cu-X distance is fixed, the Cu-Cu distance 
is variable). The upper panels show the total exchange Ji (LSDA+[7, U d = 7 eV) along with J X AFM = 4t?/[/ e ff and jf M = 
Ji — J AFM . The bottom panels show /3 4 , where /3 is the ligand's contribution to the Cu-based magnetic orbital. The WFs for 
the experimental (relaxed) geometries are shown as insets. 
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FIG. 5. Three examples of model clusters: CU3 trimer clus- 
ter as the minimal cluster for the evaluation of J\ and J 2 
(A); the pentamer cluster for calculating J 2 , with only two 
Cu 2+ and three substituted non-magnetic ions (B); and the 
tetramer cluster for calculating Ji with two magnetic Cu and 
two nonmagnetic Ms centers (C). 
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In this study, we apply the B3LYP functional on dimer 
models and vary a between 0.15 and 0.25 (a = 0.20 
corresponds to the standard B3LYP functional as im- 
plemented in Gaussian). Although we pointed out that 
dimer models are too small for calculating J\ in quanti- 
tative agreement with the periodic model, they are well 
suited for comparing the different DFT methods and pa- 
rameter sets. 66 Despite substantially different treatment 
of many-body effects in DFT+/7 and hybrid functionals, 



FIG. 6. (Color online) J 1 (N)/J 1 (N = S), J 2 (iV)/ MN = 8) 
and a(N)/a(N = 8) as a function of the chain length N. 
The bridging angle is fixed to the experimental (CuCb and 
CuBr 2 ) and optimized (CuF 2 ) values, respectively. For J 2 and 
— J2/J1, the minimal number of Cu-centers amounts to three. 
The exchange couplings are calculated with the LSDA+?7 
method with Vd — 7 eV. 
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FIG. 7. (Color online) J± and J2 of the Cu-Br clusters cal- 
culated with clusters containing two correlated and N5 un- 
corrected Mg 2+ or Zn 2+ centers. The bridging angle is fixed 
to the experimental value. The resulting exchange integrals 
are normalized to that for the Cus-octamer cluster. For the 
calculations, the LSDA+f/ method is used with Ud = 7eV. 
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FIG. 9. (Color online) Exchange integrals J\ and the ratio 
— J2/J1 of CuCb as a function of the bridging angle calculated 
with a periodic and two cluster models. The dashed vertical 
line indicates the experimental bridging angle of 93.6°. For 
the CuF2 data, see supplementary information. For the cal- 
culations, the LSDA+[7 method is used with Ud = 7eV. 
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FIG. 10. (Color online) CuBr2: the nearest-neighbor transfer 
integral t\ as a function of the bridging angle calculated with 
a periodic as well as two cluster models (CU4 and Cus). 



FIG. 8. (Color online) CuBr 2 : exchange integrals Ji and J2 
as a function of the bridging angle. A periodic as well as two 
different cluster models (CU4 and Cus) were used. The inset 
shows the ratio — J2/J1. The dashed vertical line indicates 
the experimental bridging angle of 92°. For the calculations, 
the LSDA+J7 method is used with U d = 7eV. 



produce the thermodynamical data for small molecules, 
it provides meaningful results for strongly correlated sys- 
tems such as CuX2, in line with the earlier studies. 33,67,68 
Moreover, the calculated exchange integrals are robust 
with respect to a: the exchange integrals are rather in- 
sensitive to the choice of this parameter. 



the resulting exchange integrals of all three CuX2 com- 
pounds are quite similar (Fig. 11). Thus, the B3LYP cal- 
culations confirm the LSDA+Z7 results, justify the choice 
of the free parameters in the latter approach and demon- 
strate that the unusual FM J\ coupling of CuCl2 and 
CuBr2 is not an artifact of a certain method. Despite 
the fact that B3LYP was originally constructed to re- 



V. DISCUSSION AND SUMMARY 

Our study of magnetostructural correlations in the 
CuX 2 halides reveals the crucial role of the ligand in mag- 
netic exchange. Its effect is two-fold: First, the larger 
size of CI and Br is responsible for the enhanced NNN 
coupling J2 that is assisted by the sizable overlap of lig- 
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FIG. 11. (Color online) The exchange integral J\ of the CuX 2 (X = F, CI, Br) compounds as a function of the bridging angle. 
The calculations are done for a dimer model with LSDA+?7 and Ua = 7 =b 1 eV (grey area), and with the B3LYP functional 
(a = 0.15 -0.25). 



and p orbitals along the Cu-X-X-Cu pathway. Second, 
the covalent nature of the Cu-Cl and Cu-Br bonds un- 
derlies the large ligand contribution to the magnetic or- 
bitals and, consequently, the strong FM nearest-neighbor 
(NN) coupling J\ in the broad range of bridging angles 
which could be ascribed to Hund's exchange on the ligand 
site. The tendency of covalent Cu 2+ halides to exhibit 
FM exchange along the Cu-X-Cu pathways can be illus- 
trated but also challenged by several experimental obser- 
vations. It should be emphasized that ferromagnetic NN 
coupling requires not only sizeable ferromagnetic contri- 
butions but also small transfer integrals as were found 
for CuCl2 and CuBr2. Otherwise, the AFM contribu- 
tions will outweigh the FM terms even for covalent com- 
pounds. 

Experimental data for Cu 2+ chlorides and bromides 
indeed show the robust FM NN coupling for the bridg- 
ing angles below 90°. While the < 90° regime is not 
typical for the ionic oxides and fluorides, it is abun- 
dant in covalent systems and observed, e.g., in Cu- 
based FM spin chains. 69 The FM nature of the NN cou- 
pling at = 90 — 95° is evidenced by CuCl2 and CuBr2 
themselves. 42,43,45,58 However, larger values are less 
common and require geometries other than the edge- 
sharing CuX2 chains considered in the present study. 

The angles of > 95° are only found in edge- 
sharing dimers and corner-sharing chains. Moreover, the 
respective experimental situation is rather incoherent. 
In (CuBr)LaNb 2 7 and (CuCl)LaTa 2 7 , the corner- 
sharing geometry with > 100° indeed leads to the FM 
exchange, although with a tendency towards AFM ex- 
change at > 108 - 109° (Refs. 70 and 71). By contrast, 
the CU2CI6 dimers may reveal the AFM exchange even 
at ~ 95.5°, as in LiCuCl 3 • 2H 2 (Ref. 72) or TlCuCl 3 
and KCUCI3 (Ref. 73) where the latter exhibit transfer 
integrals that are 3.5 times larger as that in CuCl2. On 
the other hand, similar CU2CI6 dimers with the same 
bridging angle of 6 ~ 95.5° in the spin-ladder compound 
IPA-CUCI3 feature the sizable FM intradimer coupling. 74 

These experimental examples show that the bridging 
angle may not be the single geometrical parameter 
determining the Cu-X-Cu superexchange. Details of 
the atomic arrangement are important even for Cu 2+ 



oxides, 20,67 whereas in more covalent systems this effect 
is likely exaggerated because interactions involve specific 
orbitals, so that each bond determines the orientation of 
other bonds around the same atom. We have pointed out, 
that such magnetostructural correlations, essential for 
understanding the magnetic behavior and for the search 
of new interesting materials, can nicely be investigated 
with cluster models. In particular in case of intricate 
crystal structures clusters enable studying effects of each 
structural parameters separately, while for periodic mod- 
els only a set of parameters can be modified at once. 

On a more general side, our results identify the Cu- 
X-Cu pathways as the leading mechanism of the short- 
range exchange in Cu 2+ halides. The fact that the mag- 
netostructural correlations weakly depend on the proce- 
dure of varying (Fig. 3) entails the minor role of direct 
Cu-Cu interactions, because the coupling always evolves 
in a similar fashion, no matter whether the Cu-Cu dis- 
tance is fixed or varied. Therefore, the nature of the 
ligand is of crucial importance, and affects the Cu-X-Cu 
hopping along with the FM contribution, presumably re- 
lated to the Hund's coupling on the ligand site. 64 In ionic 
systems, the nearest neighbor hopping increases with the 
bridging angle and dominates over the small FM contri- 
butions, thus leading to the conventional GKA behavior. 
However, the GKA behavior may be strongly altered in 
covalent compounds, as shown by our study and previ- 
ously argued in model studies on the effect of side groups 
and distortions. 20,21,67 

From the computational perspective, magnetic model- 
ing of chlorides and bromides is generally challenging. 
Although these compounds are still deep in the insu- 
lating regime, far from the Mott transition (ti <C J7 e ff, 
see Table I), the sizable hybridization of ligand states 
with correlated Cu 3d orbitals challenges the DFT+U 
approach, with correlation effects restricted to the d 
states. The microscopic evaluation of magnetic couplings 
in Cu 2+ chlorides and bromides indeed leads to large 
uncertainties. 71,75 Hybrid functionals, on the other hand, 
tend to overestimate magnetic exchange couplings 33 and 
provide a working, but empirical solution to the prob- 
lem of strongly correlated electronic systems. This calls 
for the development and application of alternative tech- 
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niques, as for instance ab initio quantum-chemical cal- 
culations, appropriately accounting for strong electron 
correlations. Since the wavefunction-based quantum- 
chemical calculations are presently restricted to finite sys- 
tems, they require the construction of appropriate clus- 
ters. This task has been successfully accomplished in our 
work. We have demonstrated that relatively small clus- 
ters with a low number of correlated centers are capable 
of reproducing the results obtained for periodic systems, 
and provide adequate estimates of the magnetic exchange 
even for the long-range Cu-X-X-Cu interactions. 

In summary, we have studied magnetostructural corre- 
lations in the family of CuX2 halides with X = F, CI, and 
Br. Our results show substantial differences between the 
ionic CuF2 and largely covalent CuCl2 and CuBr2. The 
fluoride compound behaves similar to Cu 2+ oxides, and 
shows weak FM exchange at the bridging angles close to 
= 90° along with the AFM exchange at > 100°. Go- 
ing from F to CI and Br leads to two major changes: i) the 
larger size of the ligand amplifies the AFM next-nearest- 
neighbor coupling J2; ii) the increased covalency of the 
Cu-X bonds results in the strong mixing between the 



Cu 3d and ligand p states, and enhances the FM contri- 
bution to the short-range nearest-neighbor coupling J\. 
We have constructed cluster models which, first, supplied 
an excellent description of local properties of the solids. 
Second, they turned out as highly valuable tool for inves- 
tigating magnetostructural correlations, e.g., they could 
be instrumental in the microscopic analysis of the co- 
valent Cu 2+ chlorides and bromides with interesting but 
still barely explored magnetism. Finally, they seem to be 
a viable approach to parameter-free quantum-chemical 
calculations of strongly correlated solids. 
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(Color online) Crystal structure of the CuX 2 (X = F, CI, Br) compounds. Edge-sharing CUX4 plaquettes form planar chains 
running along the b-axis. For CuF2, this is a fictitious structure which is introduced in order to investigate the effects of 
different ligand size on the intrachain magnetic couplings. Real CuF2 features a 2-dimensional structure of corner-sharing 
CuF4-plaquettes that would be inappropriate for such purposes. t\ and £2 denote nearest-neighbor and next-nearest-neighbor 
intrachain hopping, respectively, Uc is the interchain hopping which is expected to have a negligible effect on the intrachain 
couplings. 
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(Color online) Density of states (DOS) around the Fermi level of CuBr2, calculated with LDA. The gray shaded area 
corresponds to the total DOS and the red and blue lines belong to the partial Cu(3d) and Br(4p) DOS, respectively. 
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(Color online) Density of states (DOS) around the Fermi level of CuCb, calculated with LDA. The gray shaded area 
corresponds to the total DOS and the red and blue lines belong to the partial Cu(3d) and Cl(3p) DOS, respectively. 
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(Color online) Density of states (DOS) around the Fermi level of CuF2, calculated with LDA. The gray shaded area 
corresponds to the total DOS and the red and blue lines belong to the partial Cu(3<i) and F(2p) DOS, respectively. 
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(Color online) The nearest-neighbor coupling, Ji, of C11CI2 for different Cu-CI-Cu bridging angles, calculated with LSDA+£7 
where Ud is varied between 4-9 eV. A negative sign corresponds to ferromagnetic coupling. Varying Ud shifts the curves, 
however, has no qualitative effect on the magnetic coupling, thus, the ferromagnetic coupling that is obtained even at small 
and large bridging angles is not an artefact of the choice of Ud- The decreasing coupling strength for increasing Ud can be 
attributed to the reduced spatial extension of the Cu(3d) orbitals and, thus, the smaller overlap with the Cl(3p) orbitals. This 
entails smaller ligand contributions to the Wannier function which are shown to be proportional to the strength of 
ferromagnetic coupling. The vertical dotted line indicates the experimental bridging angle. 
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(Color online) Wannier functions (WFs) of CuBr2 for bridging angles of 75°, 90° and 115°. The WFs centered at the Cu-sites 
are strongly delocalized over the Br-ligands. At 90°, the small contributions of the neighboring Cu-sites vanish for symmetry 
reasons, moreover, the ligand contribution is smallest at this bridging angle. 
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(Color online) Magnetic coupling constant J\ and the ratio a — —J2/ Ji of CuF2 as a function of the bridging angle calculated 
with a solid model, Cu4-tetramer and Cus-octamer cluster models. The dashed vertical line indicates the optimized bridging 
angle of 102°. The singularity in the ratio a arises from the change of sign of J\ at about 100°. 
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(Color online) Supercells quadrupled along the b axis as used for the DFT+[7 calculations. The total exchange coupling 
constants Ji and J 2 were calculated from total energy differences of these three different arrangements of collinear spins. Red 
and blue spheres in the structures indicate Cu 2+ ions with spin up and spin down, respectively. Green balls are F, CI and Br 
ligands, respectively. 



